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SUMMARY 

/ 7 


Chemical reaction rate coefficients for the dissociation of C0 2 and CO are 
estimated hy use of collision theory and are compared with experimental results 
obtained at low temperatures. The effect of varying the number of square terms 
of energy effective in collisions on the rate coefficients is examined. The 
estimated rate coefficients are employed to study the interdependent chemical 
rate processes for the dissociation of C0 2 and its components (including five 
chemical species) coupled with the fluid flow behind a normal shock wave. 
Solutions of the differential equations are obtained for shock speeds up to 
10 km/sec at 10" 2 * to 10" 4 standard atmospheric density by three methods. Thus 
results of two simplified methods (one of which leads to a closed form solution) 
are compared with those of an exact method. Chemical relaxation effects for 
shock waves in C0 2 are compared with those in N 2 and 0 2 . Results for C0 2 are 
presented in the form of flow-field profiles of pressure, temperature, density, 
enthalpy, and species concentration, and in the form of a relaxation distance 
behind the shock as a function of shock speed and ambient density. 


INTRODUCTION 


Entry into planetary atmospheres having C0 2 components enhances the interest 
in the dissociative reactions of C0 2 . Although the estimated C0 2 concentration 
in planetary atmospheres is not large (estimates by Kaplan (ref. 2) are of the 
order of 10 percent for Venus and less than that for Mars), it is important to 
understand the chemical kinetics for C0 2 dissociation as well as for mixtures of 
C0 2 with other gases . 


1 This study has also been discussed in reference 1. The information con- 

tained in that paper has been considerably expanded in the present work. Addi- 

tional discussion of physical phenomena, methods of solution, and additional 

results are included in this paper. U ^ 
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A considerable amount of work has been done on phenomena in C0 2 in chemical 
equilibrium related to planetary atmosphere entry. For example , Yee, Bailey, and 
Woodward (ref. 3) have measured the heat transfer and have made equilibrium shock 
wave calculations in C0 2 ; Thomas (ref. t) has estimated the transport properties 
for equilibrium C0 2 ; Hoshizaki (ref. 5) has calculated the heat transfer for high 
speed flight in pure C0 2 assuming chemical equilibrium in the flow field; and 
Woodward (ref. 6 ) has calculated the equilibrium properties behind normal shocks 
in C0 2 - N 2 mixtures . 

In the present paper, attention is focused on the nonequilibrium dissoci- 
ative relaxation of C0 2 behind a normal shock wave. Our purpose is to estimate 
the nonequilibrium chemical effects on a dissociating C0 2 flow field. The 
results are of interest for C0 2 shock-tube and wind-tunnel studies and are useful 
as a first step in studying nonequilibrium effects in gas mixtures containing a 
C0 2 component. 

At present, the reaction rate coefficients for C0 2 dissociation are not well 
known. Although Brabbs , Belles, and Zlatarich (in work not yet published) and 
Gaydon and Hurle (ref. 7) have obtained some experimental results at relatively 
low temperatures, very little is known about the reactions or their rate coeffi- 
cients at the high temperatures behind strong shock waves . For that reason, the 
reactions are postulated and their rates are estimated in this paper. 


SYMBOLS 


A 

a 


a zi 


Bcp 

Bf r 


C P1 


d 


E 


c r 


Ef r 


h 


o 


i 


Avogadro’s number 

defined by equation ( 45 ) 

number of atoms of type z in species i 

coefficient in equilibrium coefficient expression ( 23 ) for reaction r 

coefficient in forward reaction rate coefficient expression ( 22 ) for 
reaction r 

specific heat at constant pressure for species i on a "per mole" basis 
average diameter of two colliding particles 

energy in equilibrium coefficient expression ( 23 ) for reaction r 

activation energy in forward reaction rate coefficient expression ( 22 ) for 
reaction r 

static enthalpy (per mass basis) 
static enthalpy (per mole basis) 

enthalpy of formation of species i (per mole basis) 
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K Cr 

k f r 

k 

M 

M* 

n i 

P 

P 

A 

R 

s 

T 

u 

X 

x 

a 

3 

e 

'H 

P 

Po 

a 


equilibrium coefficient for reaction r 
forward reaction rate coefficient for reaction r 
total number of species (eq. (7)) 

"typical 11 collision partner and molecular weight 

reduced molecular weight (l/Mj. + 1 /Mm) ~ 1 

number of moles of species i per unit mass of mixture 

transition probability or frequency factor 

static pressure 

universal gas constant 

number of classical squared terms of energy contributing to reaction 
absolute temperature 
velocity in x direction 
species 

distance behind shock 

forward stoichiometric coefficient, and exponent in equations (22) and (23) 
backward stoichiometric coefficient 

p 

— density ratio across shock wave 
Ps 

defined by equation (51) 
mass density of mixture 

standard or sea level atmosphere density (1.225x10 3 g/cm 3 ) 
constant in equation (21) 

Superscripts 

per mole basis 
equilibrium value 
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Subscripts 


a location of point in sketch (b) 

b backward 

eq equilibrium value 

f forward 

i species i; i = 1 through 5 corresponds to C0 2 , CO, 0 2 , 0, and C, 

respectively 

M pertaining to the typical collision partner 

o reference condition at sea level 

r rth reaction; r = 1 through 3 shown in equations (l8) through (20) 

s conditions immediately behind shock 

00 conditions ahead of shock 


ANALYSIS 

Differential Equations 


The flow model chosen for analysis is that of a normal shock moving at 
velocity into quiescent C0 2 at density p^- The corresponding flow field 

as seen by an observer traveling with the shock is shown in sketch (a). 



SHOCK 


Sketch (a).- Flow model. 
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The equations describing the flow field behind 
the shock are (neglecting transport phenomena) 


pu = p^ u w = constant 

(i) 

pu(du/dx) = -(dp/dx) 

( 2 ) 

u du + dh = 0 

( 3 ) 


where enthalpy of the mixture of species is 


k 

h(p, Pj n x . . . n^) = nqhp (4) 

i=i 
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and 


hi = 


T 


f* o 


Cp dT + hi 


(5) 


It is noted on the left side of equation (4) that h is assumed to be a function 
of pressure , density, and chemical composition but not of the vibrational state 
of the molecules. Thus we neglect vibrational relaxation in favor of chemical 
relaxation (the former being fast compared with the latter for C0 2 , according to 
reference 7, at temperatures up to 2700° K). The advantage of this simplifica- 
tion will be apparent subsequently when we consider specific heat. The equation 
of state of the mixture is ^ 

p = n ± (6) 

i=i 


Equations (l), (2), ( 3 ); and (4) can be combined to yield (ref. 8) 

n k 


du 

u 


dh 


p(8h/5p) i=1 dnj_ 


dn-; 


1 - u 2 


(dh/dp) - (l/p) 
-Oh/dp) ■ 


(7) 


which is readily put into the form 


du 

dx 


u 


tz 


1 - 


pu 2 1=1 


'Pi 


k 

E n i 

i= 1 H 



Pi 


f ^pi dm 
i=i fi ^ 

f s pi 

2-, n i 

i=i R 



£ Hi 


( 8 ) 


where It Is assumed for the sake of simplicity that cp^/f^ is constant at its 
fully excited classical value (5/2 for atoms, 9/2 for diatomic molecules, and 7 
for triatomic molecules (ref. 9))* 


Thus flow equations (l) through (3) are coupled to the chemical rate equa- 
tions by equation (8). The chemical rate equations for the reactions 


k 

E i x i 


1=1 


<f r k 


k b 


I>r,iXi 


’r i=i 


(9) 


can be written as (ref. 8 or 10) 


dnq 

dx 


pu 2 ^r,i - a r,i) k fj 




K, 


c r 


( 10 ) 
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If there are m types of atoms, m of equation (10) can be replaced by state- 
ments of conservation of atoms (in this case carbon and oxygen atoms) of the form 

k k 

53 a Zi n i = 53 a zi n is 

i=i i=i 

(which of course can be differentiated to give m values of (dni/dx) for use in 
eq. (8))* The boundary conditions of the flow equations (l), (2), (3)> and (8) 
and the chemical rate equation (10) are specified immediately behind the shock 
and are at x = 0, 


p 


Pfco/ £ 

(12) 

u 

= Us = 


(13) 

p 

= Ps = 

Poo^l - e) 

(lk) 

h 

= = 

ni s (c pi T s + hx°) 

(15a) 


= 

+ (l/2)(u OT 2 - u s 2 ) 

(15b) 

Hi 

c 

LI 

= (Mx ) -1 , n 2 = n 3 = n 4 = n 5 = 0 

( 16 ) 


where the strong shock approximation has been used for boundary condition (lk), 
and it has been assumed that the C0 2 does not dissociate in passing through the 
shock wave in boundary conditions (15) and (l6). 

A combination of boundary conditions (13) and (15b) with the strong shock 
approximation (ik), and the equation of state (6) yields 

6 = [2(8 Pi /r) - l]- 1 (IT) 

which is needed for the evaluation of the boundary conditions. 

Having the set of flow equations and their boundary conditions , we proceed 
to specify the chemical reactions. 

Chemical Reactions 

The following reactions are assumed for the dissociation and recombination 
of C0 2 and its constituents. 

r = 1, C0 2 + M CO + 0 + M (l8) 

r = 2, 0 2 + M 2 20 + M (19 ) 2 

r = 3, C0+M2C + 0+ M (20) 

2 The reactions t — 2 turned out to "be uniircportant in the solutions obtained* 
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The differential equations can operate on these reactions if the reaction rate 
coefficients are either known or specified. These coefficients are of the form 

kf r = 2[ PAd 2 /a(s - l):]( 2 rt^T/M*) 1 / 2 (E fr /RT) S " 1 e’^ Efr ^ RT ^ (21) 


or expressed in the Arrhenius form; 


a fr -(E f /fo) 


kf r = B fr T ” x e (22) 

The backward reaction rate coefficient is included in the equilibrium coefficient 


K Cr - (kf r /kb r ) - B Cr T e 


-(E Cr /fc) 


(23) 


This completes the set of equations used to solve the nonequilibrium flow 
behind the shock wave. It is instructive to have the equilibrium conditions 
which the nonequilibrium properties eventually approach some distance behind the 
shock wave. 


Equilibrium Conditions 


The ten equilibrium properties useful for this analysis are n, fij., n z , 113 , 
n 4 , n 5 , p, p, T> and u. They are obtained from the simultaneous solution of the 
following ten algebraic equations . By definition 


n = rii + n 2 + n 3 + n 4 + n 5 


(24) 


Statements of conservation of oxygen and carbon atoms can be expressed as (from 
eq. ( 11 )) 


nj_ +113 - n .5 + n — 2 ri] 
ni + n 2 + £5 = ni 


The equation of state is 


p = pnRT 


Equilibrium coefficients can be expressed as 

Uq ,a 

= (n 2 n 4 /ni)p = B Cl T exp(-E Cl /RT) 
Kc 2 = (n 4 /n 3 )p = B C 2 f ac 2 exp(-E C 2 /RT) 

_ a C3 a— 

Kc 3 = (h 5 h 4 /h 2 )p = B C3 T exp(-E C 3 /RT) 
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(25) 

(26) 

(27) 

(28) 

(29) 

(30) 
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and finally the flow relationships 


(31) 

(32) 


* = Pco U oo/P 
P = Poo U oo( U oo' G ) 

and 

(l/2)(u D0 2 - u 2 ) + ni hi° = £ 5 i(^Pi^ + h i°) (33) 

s i=x 


complete the set of ten simultaneous equations for the ten unknown equilibrium 
conditions. It is noted that some of the equations are nonlinear and transcen- 
dental. These ten algebraic equations were reduced to a set of three nonlinear 
transcendental equations with n 3 , n 4 , and fis as unknowns and seven equations in 
the other unknowns. The set of three equations was solved by the oldest of all 
methods, "regula falsi" (ref. 11). The other seven equations were then solved 
quite simply for the remaining seven unknowns. 


Method of Solution 


Solutions of the nonequilibrium flow field were obtained in three ways . 

These will be referred to as the complete solution, the first approximation, and 
the second approximation. The two approximate solutions were sought in an effort 
to obtain as simple a mathematical model as possible. Both approximations are 
compared with the complete solution to establish their validity. 

Complete solution.- Flow equations (l), (2), (3), and (8) were solved simul- 
taneously with the chemical reaction equations (lo) and (ll) written for the 
reactions (l8), (19), and (20) subject to boundary conditions (12) through (l6). 
The integration was performed by use of the Adams -Moulton predictor -corrector 
variable step integration scheme (ref. 12). Values for the dissociation energy 
and some of the other physical constants in the rate coefficients and equilibrium 
coefficients were either obtained or estimated from information in references 13, 
14, 15, and 16. The rate and equilibrium coefficients will be presented subse- 
quently. Computation of one example including both the nonequilibrium flow field 
and the equilibrium conditions was accomplished in approximately one minute by 
the IEM 7090 digital computer. 

Uncoupling of chemical and flow equations.- An order of magnitude analysis 
of the flow equations (see appendix) reveals that both p and h can be expected 
to be almost constant in the flow field. Examination of the flow equations (l), 
(2), and (3) shows that if any one of the variables p, h, p, and u is strictly 
constant, the other three must be constant. In what will be called the first 
approximation, we assume that this is the case. Then the flow equations can be 
omitted and the chemical processes can be studied by use of equations (10), (ll), 
and the equation of state (6). It is noted that T is still a variable because 
the n^ vary . Q 
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Closed form solution .- Another simplification can he effected by what will 
he referred to as the second approximation- In this case it is assumed that p, 
h, p, and u are constant behind the shock and that only the forward reactions of 
equations (l8) and (20) are important. Furthermore, it is assumed that the for- 
ward reaction (l8) is essentially completed before the forward reaction (20) 
begins. That is, the C0 2 is almost completely dissociated before CO begins to 
dissociate as shown in sketch (b). This assumption will be checked subsequently 



X, DISTANCE BEHIND SHOCK 

Sketch (b) 


by comparison with the results of the complete solution. With these assumptions, 
the appropriate chemical rate equation (10) and atom conservation equation (ll) 
become, for x < x a . 


dni/d x = -k fl (p/u) s n M n 1 


2112. + n 2 n.4 — 2n.n 

s 


n i + n 2 = n lg 


where 


n M = n i + n 2 + n 4 


For x > x a . 


dn 2 /dx = -k f3 (p/u) s n M n 2 


n 2 + n4 = 2n ls 


( 34 ) 

(35) 

(36) 


(37) 

(38) 

(39) 


n 2 + n 5 = n ls 


(40) 
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where 


n M = n 2 + n 4 + ns 


(4l) 


Equations (34) through (37) can he combined with the equation of state (6) and 
the forward rate coefficient equation (22) to yield a first-order differential 
equation in which n* is the only variable. 


chi}./ dx = -Bf 1 (p/pR) s af i(p/u) s n 1 (2n ls - n L ) afl exp[ - E fl (p/p) s (2n ls - nj] (42) 


with the boundary condition, at x = 0, 

n i = n i s 

If (2n ls - n 1 ) 1-a;f i in equation (42) is rewritten 


(43) 


(2n ls ) 1 " af i[l - ( ni /2n ls )] 1-a ’ :£ ’i 

and the latter egression is expanded in a power series, equation (42) can be 
integrated for af x = l/2 with the simple result for x < x a 


In 


n i 

ni q 


+ a(n x - n ls ) + 


"( n i 2 - nf s ) a3 ( n i 3 - nj s ) a m( ni m - ni s ) 


2 - 2 ! 


3-3! 


m*m; 


; an l _ e an l S 3 

2(2n ls a) + 2-4(2n ls a)‘ 


e ani (an 1 - 1) - e^Manxs - 1) 


arii 


+ 


35 . 


2*4*6(2n ls a) v 


e ani (a 2 ni 2 - 2a ni + 2) - e ani s(a 2 ni - 2an lq + 2) 


3-5 -7 

+ 2 - 4 * 6 * 8 ( 2 n x a ) 4 
s 


e an - L (a 3 n 1 3 - 3 a2 nx 2 + 6ani - 6) 


e anis (a 3 n 3 s - 3a 2 nf s + 6an ls - 6) 


3-5-7;9 

2-4-6-8- 10( 2n ls a ) 3 


e ani (a 4 n! 4 - 4a 3 nx 3 + 3 a2 nx 2 - 6an x + 6) 


an i / 4 4 
e -^(a n lc! 


. 3 3 , _ 2 2 

4a n ls + 3 a n lg 


6an lg +6) 
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2n ls P s /p\ (2nx a) 


~ Bfl J P a 8 U 


- e 


(44) 
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where 


a - - E fi (f) 8 <«) 

It is noted that x a is determined where nx/nx is some arbitrarily small 
value ( 0.001 for example). 

Equations ( 38 ) through (4l) can be manipulated in the same way to give a 
differential equation identical to equation (42) with nx replaced by n 2 , 2n ls 
replaced by 3n ls > and Bf x , and Ef x replaced by Bf 3 , ctf 3 , and Ef 3 . Its 

boundaiy condition is at x = x a 

n 2 » n lg (46). 


The solution is identical to that shown in equation (MO with the proper inter- 
change of notation; with x replaced by (x - x a ), (2n ls a) replaced by (3ni s a ) 
wherever the former appears alone in parentheses, and on the right side of equa- 
tion (44) 2ni s is replaced by 3^i s under the square root sign. 


DISCUSSION OF RESULTS 


In presenting results, we first consider reaction rate estimates. Then some 
typical solutions of the flow equations are shown, and finally relaxation dis- 
tances are presented. 


Reaction Rate Estimates 


There are few experimental results for the dissociation rates of C0 2 and CO, 
especially for the high nonequilibrium temperatures behind normal shock waves. 

For this reason, it was necessary to estimate forward reaction rates. However, 
there is no sure way to estimate dissociation rates for diatomic molecules , and 
estimates for triatomic molecules are even more uncertain. The principal uncer- 
tainties are the values of s (the number of classical squared terms of energy 
contributing to the reaction), P (the steric factor or transition probability), 
and Ef (the activation energy). The lower limit of s and the upper limit of 
P are unity. The simplest first estimate is made with these limiting values of 
s and P and with the dissociation energy for the activation energy. This approx- 
imation is equivalent to the assumption that the colliding pairs must have trans- 
lational energy greater than or equal to the dissociation energy if a reaction is 
to occur, and that a collision of any such pair always results in a reaction. 

The corresponding rate coefficient estimates for CO 2 and CO dissociation are shown 
as solid lines in figure 1 . 

Although, in fact, s is greater than unity and P is less than unity, it 
can be shown by the following discussion that for reasonable values of each, the 
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effects tend to cancel one another and the simple estimate used here is plausible. 
The upper limit of s for the collision of a C0 2 molecule with an atom is 6. 

The limit is higher for collision with another molecule. Similarly, the upper 
limit of s for collision of CO with an atom is either 2 or 3, and for collision 
with a diatomic molecule is either 3 or 5* If we assume for illustrative purposes 

that s = 5 for C0 2 dissociation and 3 for CO dissociation, while P is unity, 

the reaction rate coefficients for these two reactions are shown by dashed lines 
in figure 1. At temperatures below 20,000° K, they are higher than the corre- 
sponding solid lines by a factor of the order of 10 to 10 3 . However, if P were 
of the order of 10 " 1 to 10 " 3 , the effect would be canceled and we would come 
back to about the same estimate obtained for s = P = 1. Examination of rate 
data for oxygen dissociation Indicates that P is probably of the order of 10~ 2 
for that reaction. Thus we conclude that the simple first estimate using 
P = s = 1 leads to plausible rate coefficients in the temperature range up to 
20,000° K. At temperatures greater than 20,000° K, the rates for s > 1 and 
P = 1 differ from the estimates for s = P = 1 by less than a factor of 10. 

Thus since s is probably greater than 1, while P Is much less than 1, the 

rates predicted by setting s = P = 1 are expected to be high. Because of the 
preceding argument, it is speculated that relaxation phenomena may be even slower 
and relaxation distances may be even larger than the results of the present 
estimates will show. 

Experimental results for the reaction 

C0 2 + A -* products (47) 

have been obtained (by Brabbs, Belles, and Zlatarich) at temperatures to about 
3000° K, from which the reaction rate coefficient 

kf = 3x10 T exp ( -86/RT)cm 3 mole“ 1 sec” 1 (48) 

was obtained, which is also shown in figure 1. It is noted that the curve 
corresponding to equation (48) is in fair agreement with the present estimate for 
C0 2 dissociation in the temperature range 5^000° K to 10,000° K. 

The curve for oxygen dissociation shown in figure 1 was obtained from refer- 
ence 15 (in which the results of ref. 1 6 were extrapolated to temperatures in 
excess of 20,000° K). 

Equilibrium coefficients for the three reactions considered are shown in 
figure 2. These were obtained by an empirical fit of equation (23) to the 
numerical results of references 12 and 13* 

A summary of the constants used in the forward reaction rate coefficients 
and equilibrium coefficients Is shown In the following table. 






Reaction 

summary 




Bf r , 


Ef r , 

®c r J 


E Cr > 

r 

cm 3 

af 

± T 

k cal 

g mol 

a c r 

k cal 


g mol sec 0 K a ^ r 


ghnol 

cm 3 °K ac r 


g mol 

1 

6. 955X10 12 

0.5 

125- 7 497 

0.5 05 96x10 14 

-2.9713 

138.7454 

2 

11.9X10 20 

-1.5 

117 . 9982 

•33399X10 4 

-.6038 

119.4466 

3 

7.7206xl0 12 

• 5 

256.1742 

. 34113x10 s 

-.2366 

257-4958 


Solutions 


Before a few typical solutions are examined, it is instructive to examine 
the equilibrium values which the flow-field properties approach asymptotically 
some distance behind the shock. Both static enthalpy and pressure are almost 
constant behind the shock (see appendix). Thus we can expect that p/p s and 
h/h s are approximately unity. Indeed, for flight speeds between 5 and 10 km/sec 
in C0 2 at ambient density between 10" 2 and 10” 4 , standard atmospheric density 
p/p s and h/h s differed from unity by less than 3 percent. First, estimates of 
p and h can thus be obtained from equations (14) and (15b) directly. However, 
p and T, shown in figures 3 and 4, are quite different from and T g . It is 
seen that p is not a strong function of shock speed. Moreover, p^/p varies 
between the narrow limits 0.0472 to 0.0628 for the entire range of conditions men 
tioned previously. Equilibrium temperatures shown in figure 4 are fairly strong 
functions of both shock speed and ambient density. 

Turning now to the solutions of the chemical rate and flow equations, we 
examine first a case of a shock moving at fairly high speed (9 km/sec) into C0 2 
at low density (p^/pQ - 10“ 4 ). Flow-field quantities are shown in figure 5(a) 
and chemical species concentrations are shown in 5(b). Note that quite large 
variations in the flow-field quantities u, p, and T occur because of the 
chemical relaxation processes. 

It is interesting to observe that behind the shock, the flow accelerates and 
then decelerates because of the chemical relaxation processes. Correspondingly, 
density diminishes and then increases so that the product of density and velocity 
remains constant (eq. (l)). Generally, at lower shock speeds in air, density 
increases monotonically in the dissociative relaxation region behind the shock 
(ref. 17). For that reason, the negative density gradient behind the shock in 
C0 2 may be a little surprising. 

The conditions for which the density gradient behind the shock is positive 
or negative can be derived as follows. Just behind the shock, n 2 = n 3 = n 4 = ns =0 
and dna/dx = dns/dx = 0 (because there is no atomic oxygen to recombine and 
there is no carbon monoxide to dissociate). Using this information with the equa 
tion of state (6), boundary conditions ( 12) , (13), (l4), equations (17), (ll), and 
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(8) leads to 


du\ 

dx Js 



A A 

C P1 C P 
-2 + 


A 

R 


A 

A ^ A 

R R 


4 - hi° + h 2 ° + h 4 ° 

+ V 2 2(8p/ft - 1) 


n is [2(S Pi /R) -1]' 


( 49 ) 


As long as C0 2 is dissociating behind the shock, (dni/dx) s is negative. Then 
(du/dx) s is + or (dp/dx) s is ± if the traces in equation (49) are ±. Rearrang- 
ing the quantities in the braces and substituting numerical values for h^ 0 and 
ni, leads to the criteria 


d£\ 

dx/s 


£ 0 if u^ 2 $ 



cm 2 /sec 2 


(50) 


Applying equation (50) to the case where all species are fully excited leads to 
the approximate result (dp/dx) s > 0 if Uoo ^ 5 km/sec. For no vibrational 
excitation, the 5 is changed to about 9 km/sec. 

In figure 5> the solid lines correspond to the complete solution of the 
coupled chemical and flow equations . The long dashed lines correspond to the 
first approximation in which the chemistry and flow are uncoupled and only the 
chemical rate equations are solved. The short dashed lines correspond to the 
closed -form solution of the second approximation in which only the forward reac- 
tions are considered. For this flight condition, neither approximation repre- 
sents the temperature or concentration profiles of the complete solution very 
well. Evidently the coupling between the flow field and the chemical rate proc- 
esses is important; and, as can be seen in figure 5(b), C0 2 is not completely 
dissociated before CO begins to dissociate, contrary to the assuBption made in 
the second approximation. It is particularly noteworthy in figure 5(a) that the 
flow field is still far from equilibrium 10 centimeters behind the shock. 

It is interesting to observe the effect of a large change in reaction rate. 
In figures 6(a) and 6(b), the solid lines correspond to the complete solution 
shown previously in figure 5 • The dashed lines are the result of changing s in 
reactions (l8) and (20) to 5 and 3, respectively, leaving P fixed at unity. 

This corresponds loosely to increasing the reaction rates by about two or three 
orders of magnitude. It is seen that equilibrium is approached more rapidly for 
the latter condition (although relaxation distance is only changed about one 
order of magnitude) as would be expected. However, for the reasons mentioned 
previously, the solid lines are believed to be a more realistic representation of 
the flow field. 
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Figures 7(a) and 7(b) show flow -field and concentration profiles for a lower 
flight speed at the same ambient density. Here the first and second approxima- 
tions give almost identical results and both are a reasonable representation of 
the temperature and concentration profiles of the complete solution. The results 
indicate that only reaction (l8) was of any consequence at this speed. 

A result of high-speed flight (u^ = 10 km/sec) at one hundred times higher 
density (p^Po = 10 -2 ) is shown in figures 8(a) and 8(b). Again for high speed, 
neither approximation represents the complete solution very well. Here it is 
noted that most of the chemical relaxation has taken place in a distance of 
0.1 cm behind the shock. 

In figures 9(a) and 9(b), the flow field is that of a low shock speed at a 
high ambient density. Here again, the dissociation of C0 2 was the only signifi- 
cant reaction, as may be expected. 3 Both approximations were reasonable 
representations of the flow field. 


Comparison of C0 2 With Other Gases 


For strong shock waves, the flow behind the shock can be characterized by 
two free -stream variables, velocity and ambient density. For a fixed free -stream 
velocity and ambient density, we will compare results in C0 2 with those in air, 
nitrogen, and oxygen. 

Insofar as equilibrium conditions behind shock waves are concerned, the 
pressure is very nearly the same for all of these gases. However, considerably 
lower equilibrium temperature, density, and mole fraction of ions or electrons 
exist behind high-speed shock waves in C0 2 than in air. 

With regard to nonequilibrium processes, it is pointed out that in terms of 
dissociation energy, C0 2 is more stable than 0 2 but less stable than W 2 • Also, 

CO is more stable than N 2 . Thus for high-speed shock waves for which both C0 2 
and the resulting CO dissociate, we would expect relaxation to proceed more 
slowly in C0 2 than in 0 2 , N 2 , or air. Conversely, for low -speed shock waves 
behind which C0 2 dissociates but CO does not, we would expect relaxation to 
proceed more rapidly in C0 2 than in N 2 . 

The nonequilibrium flow -field profiles behind shocks moving at 9 km/sec at 
ambient density P^/po = 10“ 4 in pure C0 2 and pure N 2 are compared in figure 10. 
Reaction rates for N 2 dissociation were obtained from reference l8. Relaxation 
in C0 2 is slower than in N 2 , as was expected. The initial dip in the density 
profile appears to be missing in the N 2 flow. It is actually there, however, but 
is small and almost at the origin. 

A comparison of C0 2 with W 2 and 0 2 can be summarized in terms of the degree 
of completion q (which is defined subsequently) of the chemical relaxation. 

That is, at a distance of 10 centimeters behind the shock wave (again 

3 In figure 9(b), the concentration of 0 is slightly less than that of CO 
because some of the atomic oxygen recombined to form a trace of 0 2 . 
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Uqq = 9 km/sec , P^/pQ = 10” 4 ), the chemical relaxation is about 60 percent com- 
plete in C0 2 , but is about 85 percent complete in N 2 , -while at a distance of only 
0.05 cm behind the shock in 0 2 , the chemical relaxation is 99 percent complete. 

The above results will be summarized in a general way subsequently. Before 
drawing our conclusions, however, it is worthwhile to examine the results of 
numerous solutions from the standpoint of relaxation distances. 


Relaxation Distances 


For the present purposes, the relaxation distance is defined as the distance 
behind the shock at which the change in number of moles that has occurred is a 
specified fraction of the total change necessary to achieve equilibrium. That is, 
the relaxation distance is that distance behind the shock at which 

k 

E n i " n is 

11 = FT (51) 

_ ni s 

i=i 

where rj is a specified value less than unity. Thus rj is an index of how far 
the chemical reactions have gone toward completion. Figures 11 through Ik are 
contour charts of relaxation distance as a function of shock speed and ambient 
density for r\ = 0.1, 0.5, 0.8, and 0.95; respectively. It can be seen that the 
relaxation distance for 10 -percent reaction completion is small, varying from 
10~ 5 to 10“ 2 cm as p /p~ varies from 10“ x to 10 -4 . On the other hand, relaxa- 
tion distances for 95 -percent completion can be large, varying from 10 to 10 cm 
as P^/pq varies from lO" 1 to 10” 4 . 

CONCLUDING REMARKS 


A study was made of the dissociative relaxation of carbon dioxide and its 
components behind normal shock waves moving at speeds between 5 and 10 km/sec at 
ambient densities from 10~ 2 to 10" 4 times standard atmospheric density. Because 
of a lack of high -temperature chemical reaction-rate data, it was necessary to 
use collision theory for estimating rate processes. The analysis was simplified 
somewhat and it was assumed that vibrational relaxation could be neglected. It 
was also assumed that the specific heat of a pure species was constant at its 
classical fully excited value. The latter assumption is quite reasonable for 
single particle species C and 0 in the regime of the analysis, but less valid for 
the species CO and C0 2 * Certainly when more is known of reaction rates, It would 
be appropriate to refine the other assumptions used in the study. 

The dissociation and recombination reactions for C0 2 , 0 2 , and CO were 
considered, in the analysis, and the results indicated that the dissociation and 
recombination of O2 could have been neglected. The analysis results also showed 
that at large distances behind the shock for shock speeds of about 6 km/sec, C0 2 
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was almost completely dissociated but CO did not dissociate. On the other hand, 
at shock speeds of 10 km/sec almost everything was dissociated and the end prod- 
ucts consisted principally of atomic carbon and oxygen. It should be noted, how- 
ever, that ionizing reactions were not considered in the analysis and that, at 
10 km/sec, some ionization may be present. 

Strong coupling was indicated between the chemical relaxation and the flow 
field in the nonequilibrium part of the flow. The coupling effects were strong 
for high-speed shocks regardless of density level. The coupling for high-speed 
shocks causes density to first decrease and then increase toward the equilibrium 
value. Particle velocity behaves in the reverse way. A simple criterion was 
developed to show conditions for which the density gradient behind the shock wave 
is positive or negative. 

Two approximations were employed to obtain solutions. The first of these 
uncoupled the chemical effects from the flow field. The second, which uncoupled 
the chemical effects and in addition considered only the dissociation reactions, 
led to a closed-form solution. Both approximations take into account the varying 
temperature in the relaxation region. Comparison of the results with those of 
the complete solutions indicated that both approximations gave about the same 
result for temperature and concentration profiles at shock speeds up to 6 km/sec 
and both were a reasonable representation of the complete solution. At shock 
speeds of 10 km/sec, however, neither approximation represented the complete 
solution very well. 

Results show that flow fields behind high-speed shock waves are farther from 
chemical equilibrium in C0 2 than in N 2 or 0 2 . 

Finally, relaxation distances behind the shock were large for low density. 
For ambient density of the order of 10 -4 standard density, the relaxation 
distance for reactions to go 95 percent toward completion was of the order of 
10 meters . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., March 19, 19^3 
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APPENDIX 


ORDER OF MAGNITUDE DISCUSSION 


In the strict sense, all of the flow-field quantities p, u, p, h, T and 
species concentration vary in the reacting flow field behind the shock. Some of 
these quantities, however, will vary more than others. The purpose of the follow- 
ing discussion is to show the relative variation of the quantities p, u, p, and 
h in the flow field behind the shock. 

At some distance x behind the shock, the ratio of local density to the 
density immediately behind the shock is defined as g. 


(Al) 


Thus the change in p from p g is 


P " Ps AP , 

= — = g - ! 


p s p s 


The corresponding change in u is from equation (l) 


Au 1 . 

u s _ g " 


(A2) 


(A3) 


Similarly, from equation (2) the change in p (noting that p s ~ p^u^ 2 = p s u s u co ) 
is 


AP _ e 
Ps 



(A*0 


where e = p /p s . Finally, the change in h if we assume that h^ « h g and 
Uoo 2 » u s and make use of the integrated form of equation (3) is approximately 


Ah 

hs 




(A5) 


Assuming for exemplary purposes that g may be as high as 2, and that 
e x 10 -1 , the corresponding changes in density, velocity, pressure, and enthalpy 
(from eqs . (A2) through (A5)) are 100, 50, 5, and 0-75 percent, respectively. 
Similarly, if g is as small as l/2, we arrive at the same qualitative result. 
Thus we conclude that pressure and enthalpy do not vary greatly in the chemical 
relaxation region, although density and particle velocity vary a great deal. 
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C0 + 0,(S = 5, P = l) 



Figure 1. - Forward reaction rate coefficients 
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Figure 3.- Equilibrium values of density 'behind normal shock. 
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X, DISTANCE BEHIND SHOCK, cm 

(a) Velocity and thermodynamic properties. 

Figure 5 • _ Nonequilibrium flow-field profiles (u^ = 9 km/sec, 0„/p O - 10-4 )• 
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DISTANCE BEHIND SHOCK, cm 

(t) Species concentration. 

Figure 5.- Concluded. 
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DISTANCE BEHIND SHOCK, Cm 
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X, DISTANCE BEHIND SHOCK, cm 

(a) Velocity and. thermodynamic properties. 

Figure 7 .- Nonequilibrium. flow-field profiles (u^ = 6 km/sec, P^/pQ = 10 
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(t) Species concentration. 
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COMPLETE SOLUTION 
FIRST APPROXIMATION 
SECOND APPROXIMATION 
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(t>) Species concentration. 
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X, DISTANCE BEHIND SHOCK, cm 

(a) Velocity and thermodynamic properties. 

Figure 9- “ Nonequilfbrium flow-field profiles (u w = 6 km/sec, P„/»o ■ 1Q - 2 )- 
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Figure 9- - Concluded. 
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DISTANCE BEHIND SHOCK (cm) 

Comparison between C0 2 and N 2 = 9 km/sec, P^/pQ = 10 -4 ). 








